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Dihadron angular correlations in d+Au collisions at = 200 GeV are reported as a function 

of the measured zero-degree calorimeter neutral energy and the forward charged hadron multiplicity 
in the Au-beam direction. A finite correlated yield is observed at large relative pseudorapidity (A 77 ) 
on the near side (i.e. relative azimuth Acj ~ 0). This correlated yield as a function of Arj appears to 
scale with the dominant, primarily jet-related, away-side (A0 ~ tt) yield. The Fourier coefficients of 
the Acp correlation, Vn = (cos nA(j >), have a strong Ar; dependence. In addition, it is found that V\ is 
approximately inversely proportional to the mid-rapidity event multiplicity, while V 2 is independent 
of it with similar magnitude in the forward (d-going) and backward (Au-going) directions. 

PACS numbers: 25.75.-q, 25.75.Dw 


Relativistic heavy-ion collisions are used to study 
quantum chromodynamics (QCD) at high energy den¬ 


sities at the Relativistic Heavy Ion Collider (RHIC) and 
the Large Hadron Collider (LHC) [TH5]. Final-state par- 
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tide emission in such collisions is anisotropic, quantita¬ 
tively consistent with hydrodynamic flow resulting from 
the initial-state overlap geometry El E]. Two-particle 
correlations are widely used to measure anisotropic flow 
and jet-like correlations [5]. A near-side long-range cor¬ 
relation (at small relative azimuth Acj) and large rela¬ 
tive pseudorapidity Arj), called the “ridge,” has been ob¬ 
served after elliptic flow subtraction in central heavy-ion 
collisions at RHIC and the LHC [9HI1]. It is attributed 
primarily to triangular flow, resulting from a hydrody¬ 
namic response to initial geometry fluctuations usms]. 

As reference, p+p, p+A and d-|-Au collisions are of¬ 
ten used to compare with heavy-ion collisions. Hydro¬ 
dynamics is not expected to describe these small-system 
collisions. However, a large Arj ridge has been observed 
in high-multiplicity p+p m and p-fPb [MS] colli¬ 
sions at the LHC after a uniform background subtrac¬ 
tion. The similarity to the heavy-ion ridge is suggestive 
of a hydrodynamic description of its origin, in conflict 
with early expectations. Indeed, hydrodynamic calcula¬ 
tions with event-by-event fluctuations can describe the 
observed ridge and attribute it to elliptic flow EH ESj ■ 
Other physics mechanisms are also possible, such as the 
color glass condensate where the two-gluon density is en¬ 
hanced at small Acj) over a wide range of Arj pH - ESj . or 
quantum initial anisotropy m- 

Furthermore, a back-to-back ridge is revealed by sub¬ 
tracting dihadron correlations in low-multiplicity p-fPb 
from those in high-multiplicity collisions at the LHC [T9|- 
A similar double ridge is observed in d-|-Au collisions 
at RHIC by PHENIX within 0.48 < \Ari\ < 0.70 us¬ 
ing the same subtraction technique EH]. A recent STAR 
analysis has challenged the assumption of this subtrac¬ 
tion procedure that jet-like correlations are equal in high- 
and low-multiplicity events |29] . It was shown that the 
double ridge at these small-to-moderate Arj has a sig- 
niflcant contribution from residual jet-like correlations 
despite performing event selections via forward multi¬ 
plicities ES]- A recent PHENIX study of large Ap cor¬ 
relations, without relying on the subtraction technique, 
suggests a long-range correlation consistent with hydro- 
dynamic anisotropic flow [30j . In order to further un¬ 
derstand the underlying physics mechanism, here in this 
Letter, we present our results on long-range (large Ap) 
correlations in d-|-Au collisions at = 200 GeV as a 

function of A 77 and the event multiplicity. The large ac¬ 
ceptance of the STAR detector is particularly well suited 
for such an analysis over a wider range in Ap. 

The data were taken during the d-|-Au run in 2003 by 
the STAR experiment EH EH • The details of the STAR 
detector can be found in Ref. ESj- Minimum-bias d-l-Au 
events were triggered by coincidence of signals from the 
Zero Degree Calorimeters (ZDC) [33] the Beam- 
Beam Counters (BBC) EH- Particle tracks were recon¬ 
structed in the Time Projection Chamber (TPC) ES] and 
the forward TPC (FTPC) ESj- The primary vertex was 


determined from reconstructed tracks. In this analysis, 
events were required to have a primary vertex position 
l-Zvtxl < 50 cm from the TPC center along the beam axis. 
TPC(FTPC) tracks were required to have at least 25(5) 
out of the maximum possible 45(10) hits and a distance 
of closest approach to the primary vertex within 3 cm. 

Three measurements were used to select d-|-Au events: 
neutral energy by the ZDC and charged particle multi¬ 
plicity within —3.8 < p < —2.8 by the FTPC EH ISHj 
both in the Au-beam direction, and charged particle 
multiplicity within \p\ < 1 by TPC. Weak but posi¬ 
tive correlations were observed between these measure¬ 
ments; the same event fraction defined by these measures 
corresponded to significantly different d-|-Au event sam¬ 
ples ESI- In this work we study 0 - 20 % high-activity and 
40-100% low-activity collisions according to each mea¬ 
sure. 

The pairs of particles used in dihadron correlations are 
customarily called the trigger and the associated particle. 
Two sets of dihadron correlations are analyzed: TPC- 
TPC correlations where both the trigger and associated 
particles are from the TPC (|? 7 | < 1), and TPC-FTPC 
correlations where the trigger particle is from the TPC 
but the associated particle is from either the FTPC-Au 
(-3.8 < p < -2.8) or FTPC-d (2.8 < p < 3.8). The 
Pt ranges of the trigger and associated particles are both 
1 < Pt < 3 GeV/c. The associated particle yields are 
normalized per trigger particle. The yields are corrected 
for the TPC and FTPC associated particle tracking ef¬ 
ficiencies of 85% ± 5% (syst.) and 70% ± 5% (syst.), re¬ 
spectively, which do not depend on the event activity in 
d-l-Au collisions EHI32]. 

The detector non-uniformity in Acj) is corrected by the 
event-mixing technique, where a trigger particle from one 
event is paired with associated particles from another 
event. The mixed events are required to be within 1 cm in 
•Zvtx, with the same multiplicity (by FTPC-Au or TPC) 
or similar energy (by ZDC-Au). The mixed-event corre¬ 
lations are normalized to 100% at Ap = 0 for TPC, and 
at ±3.3 for FTPC-d and FTPC-Au associated particles, 
respectively. 

Two analysis approaches are taken. One is to analyze 
the correlated yields after subtracting a uniform combi¬ 
natorial background. The background normalization is 
estimated by the Zero-Yield-At-Minimum (ZYAM) as¬ 
sumption P ETj- ZYAM is taken as the lowest yield 
averaged over a Acj) window of 7 r /8 radian width, after 
the correlated yield distribution is folded into the range 
of 0 < Acj) < TT. The ZYAM systematic uncertainty is 
estimated by the yields averaged over windows of half 
and three half the width. We also fit the Acj) correla¬ 
tions by two Gaussians (with centroids fixed at 0 and tt) 
plus a pedestal. The fitted pedestal is consistent with 
ZYAM within the statistical and systematic errors be¬ 
cause the near- and away-side peaks are well separated 
in d±Au collisions. The systematic uncertainties on the 
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correlated yields are taken as the quadratic sum of the 
ZYAM and tracking efficiency systematic uncertainties. 
The other approach is to analyze the Fourier coefficients 
of the A(j> correlation functions, Vn = {cosnAc/)). No 
background subtraction is required. Systematic uncer¬ 
tainties on the Fourier coefficients are estimated, by vary¬ 
ing analysis cuts, to be less than 10 % for Vi and V 2 , and 
smaller than the statistical errors for V 3 . 

Figure [2 shows the ZYAM-subtracted correlated yields 
as a function of in ZDC-Au low- and high-activity 
d-|-Au collisions. The TPC-TPC correlation at large Arj 
is shown in panel (a), whereas the TPC-FTPC correla¬ 
tions are shown in panels (b) and (c) for Au- and d-going 
directions, respectively. The ZYAM statistical error is 
included as part of the systematic uncertainty drawn in 
Fig.[2because it is common to all Acj) bins. No difference 
is observed in TPC-TPC correlations between positive 
and negative A 77 , so they are combined in Figj^a). The 
away-side correlated yields are found to be larger in high- 
than low-activity d-|-Au collisions for TPC and FTPC-Au 
correlations. The opposite behavior is observed for the 
FTPC-d correlations, Fig. [^c). 

On the near side, the correlated yields are consistent 
with zero in the low-activity events and, in FTPC-d, in 
the high-activity events as well. (Note that the yield 
value cannot be negative because of the ZYAM assump¬ 
tion.) In contrast, in TPC and FTPC-Au, finite corre¬ 
lated yields are observed in high-activity events. A sim¬ 
ilar result was observed by PHENIX [3^. In Fig. the 
event activity is determined by ZDC-Au. For event ac¬ 
tivity determined by FTPC-Au or TPC multiplicity, the 
data are qualitatively similar. In Table the correlated 
yields integrated over the near side (|A^| < tt/S) and 
the away side (|A(/) — 7 r| < tt/S), normalized by the in¬ 
tegration range, are tabulated together with the ZYAM 
magnitude for low- and high-activity events determined 
by the various measures. 

For trigger particles in our px range of 1 < pt < 
3 GeV/c, the away-side correlation in d-l-Au collisions is 
expected to be dominated by jet-like correlations [35] ■ In¬ 
specting the near-side correlation amplitude at large Ap, 
any possible non-jet, e.g. anisotropic flow, contributions 
on the away side should be order of magnitude smaller. 
Perhaps the observed away-side dependence on ZDC-Au 
event activity arises from a correlation between jet pro¬ 
duction and the forward beam remnants. Or, the under¬ 
lying physics may be more complex; for example the op¬ 
posite away-side trends in the Au- and d-going directions 
may arise from different underlying parton distributions 
in high- and low-activity collisions. The finite correlated 
yield on the near side is, on the other hand, rather sur¬ 
prising because jet-like contributions should be minimal 
at these large Arj distances. Hijing simulation [38] of 
d-|-Au collisions indicates that jet correlations within our 
Pt range after ZYAM background subtraction is consis¬ 
tent with zero at jApj >1.5. 


To study the Ap dependence of the correlated yields 
in the TPC and FTPC, the correlation data are divided 
into multiple Arj bins. In Fig. [^a), the near- and away- 
side correlated yields are shown as a function of A 77 . To 
avoid auto-correlations, we have used ZDC-Au for event 
selections for both the TPC and FTPC correlation data. 
Unlike in Fig. the ZYAM statistical errors are depen¬ 
dent of Ar] and are therefore included in the statistical 
error bars of the data points. The away-side correla¬ 
tion shape, noticeably concaved for TPC, is presumably 
determined by the underlying parton-parton scattering 
kinematics. On the near side, finite correlated yields are 
observed at large Ar] on the Au-going side in all bins, 
while the yields are consistent with zero on the d-going 
side. As aforementioned, similar results have been previ¬ 
ously observed in heavy-ion mn], p+p [m, and p-fPb 
collisions There, the trigger and associated par¬ 

ticles were taken from the same rj region. As a result, the 
correlated yields were approximately uniform in Ar] [39j , 
and were dubbed the “ridge.” In the three groups of cor¬ 
relation data in Fig.|^a), the trigger particles come from 
the TPC, but the associated particles come from different 
r] regions. Significant differences in pair kinematics result 
in the steps at Ar] = ±2 even though their Ar] gaps are 
similar. Despite this, for simplicity, we refer to the large 
Ar] correlated yields in our data also as the “ridge.” 

In order to elucidate the formation mechanism of the 
ridge, we study in Fig. Sb) the ratio of the near- to 
away-side correlated yields. Because the ZYAM value is 
common for the near and away side, its statistical error 
is included as part of the systematic uncertainty; this 
part of the systematic uncertainty is uncorrelated be¬ 
tween Ar] bins. While the large peak at Arj ^ 0 is due 
to the near-side jet, the ratio at Ar] < —1 is rather in¬ 
sensitive to Ar], whether the correlations are from TPC 
or FTPC-Au. A linear fit (dashed-line in Fig. [^b)) to 
those data points at Ar] < — 1 yields a slope parameter 
of —0.023 ± O.OlOlg g^Q with y^/ndf = 2.6/3, indicating 
that the ratio is consistent with a constant within one 
standard deviation. The rather constant ratio is remark¬ 
able, given the nearly order of magnitude difference in 
the away-side jet-like correlated yields across At] = —2 
due to the vastly different pair kinematics. Since the 
away-side correlated yields are dominated by jets [55] . 
the finite, A 77 - independent ratio at Aij < — 1 may sug¬ 
gest a connection between the near-side ridge and jet 
production, even though any possible jet contribution to 
the near-side ridge at jAyyl > 1 should be minimal. On 
the other hand, the near-side ridge does not seem to scale 
with the ZYAM value, which represents the underlying 
background. A linear fit to the ratio of the near-side cor¬ 
related yield over ZYAM in the same A 77 < — 1 region 
gives a slope parameter of 6.5 ± I. 6 / 2 ’} x 10“^, signifi¬ 
cantly deviating from zero. 

The correlated yields discussed above are subject to the 
ZYAM background subtraction. Another way to quantify 
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FIG. 1: Correlated dihadron yield, per radian per unit of pseudorapidity, as a function of Acp in three ranges of At; in d+Au 
collisions. Shown are both low and high ZDC-Au activity data. Both the trigger and associated particles have 1 < Pt < 3 GeV/c. 
The arrows indicate ZYAM normalization positions. The error bars are statistical and histograms indicate the systematic 
uncertainties. 


TABLE I: Near- (|A(^| < tt/S) and away-side (|A(^ —7r| < tt/S) correlated yields and ZYAM background magnitude, per radian 
per unit of pseudorapidity, at large A?; in low- and high-activity d-fAu collisions. Positive(negative) rj corresponds to d(Au)- 
going direction. Both the trigger and associated particles have 1 < pt < 3 GeV/c. All numbers have been multiplied by 10 ^. 
Errors are statistical except the second error of each ZYAM value which is systematic and applies also to the corresponding 
near- and away-side yields. An additional 5% efficiency uncertainty applies. 


Event 

activity 

Event 1.2 < |Ap| < 1.8 

selection ZYAM near away 

Event —4.5 < Ap < —2 2 < Ap < 4.5 

selection ZYAM near away ZYAM near away 

40-100% 

0 -20% 

ZDG 1896±7tj3 346±5 

3043±11/26 53±7 456±7 

ZDG 978±2ti 2±1 55±1 361±lt2 1±1 38±1 

m&±4+\ 10±2 70±2 438±2t^ 1±1 31±1 

40-100% 

0 -20% 

FTPG 1324±7/6 7±4 347±5 

3468±10t5 43±6 429±7 

TPC 636±2ti 6±1 59±1 309±2tj 3±1 45±1 

1899±3/5 15±2 75±2 445±lt3 2±1 27±1 


the ridge is via Fourier coefRcients of the azimuthal cor¬ 
relation functions without background subtraction. Fig¬ 
ure]^ shows the second harmonic Fourier coefficient (V 2 ) 
as a function of Arj for both high and low ZDC-Au energy 
collisions. The V 2 values are approximately the same in 
high- and low-activity collisions at large Ap. Both de¬ 
crease with increasing jApj from the small Ap, jet dom¬ 
inated, region to the large Ap, ridge, region by nearly 
one order of magnitude. The Ap behavior of V 2 , a mea¬ 
sure of modulation relative to the average, is qualitatively 
consistent with the Ayy-dependent ratio of the near-side 
correlated yield over ZYAM. One motivation to analyze 
correlation data using Fourier coefficients is their inde¬ 
pendence of a ZYAM subtraction procedure. One way for 
V 2 to develop is through final-state interactions which, if 
prevalent enough, may be described in terms of hydro- 
dynamic flow. If V 2 is strictly of a hydrodynamic elliptic 
flow origin, the data would imply a decreasing collective 
effect at backward/forward rapidities that is somehow 
independent of the activity level of the events. 

To gain further insights, the multiplicity dependencies 
of the first, second and third Fourier coefficients bi, V 2 
and V 3 are shown in Fig. Three Ap ranges are pre¬ 
sented for FTPC-Au, TPC, and FTPC-d correlations, re¬ 


spectively. Results by both the ZDC-Au and FTPC-Au 
event selections are shown, plotted as a function of the 
corresponding measured charged particle pseudorapidity 
density at mid-rapidity dNch/dij. The absolute value of 
the Vi parameter in each Ap range varies approximately 
as {dNch/dri)~^ (see the superimposed curves). This is 
consistent with jet contributions and/or global statistical 
momentum conservation. On the other hand, the V 2 pa¬ 
rameter in each Ap range is approximately independent 
of dNch/dr] over the entire measured range (the dashed 
lines are to guide the eye). Similar behavior of V 2 is 
also observed in p-|-Pb collisions at the LHC [I3IM1I1T]. 
Figure 1^ shows that the V 3 values are small and mostly 
consistent with zero, except for TPC-TPC correlation at 
the lowest multiplicity. 

In d-|-Au collisions, dihadron correlations are domi¬ 
nated by jets, even at large Ap, where the away-side 
jet contributes [35] • The behavior of Vi suggests that 
the jet contribution to V„ is diluted by the multiplicity. 
The similar V 2 values and Ap dependencies in different 
multiplicity collisions are, therefore, rather surprising. 
In order to accommodate a hydrodynamic contribution, 
there must be a coincidental compensation of the reduced 
jet contribution with increasing multiplicity, over the en- 
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FIG. 2: The Arj dependence of (a) the near- (|A(^| < tt/S) 
and away-side (|A</) — 7r| < rr/S) correlated yields, and (b) 
the ratio of the near- to away-side correlated yields in d-|-Au 
collisions. Positive (negative) rj corresponds to d(Au)-going 
direction. Only high ZDC-Au activity data are shown. The 
error bars are statistical and histograms indicate the system¬ 
atic uncertainties (for Ar; > 2 in (b) the lower bound falls 
outside the plot). The dashed curve in (b) is a linear fit to 
the A ?7 < — 1 data points. 



Ar| 


FIG. 3: The At; dependence of the second harmonic Fourier 
coefficient, V 2 , in low and high ZDG-Au activity d-|-Au colli¬ 
sions. The error bars are statistical. Systematic uncertainties 
are 10% and are shown by the histograms, for clarity, only for 
the high-activity data. 


tire measured multiplicity range, by an emerging, non-jet 
contribution, such as elliptic flow. 

Whether or not a finite correlated yield appears on 
the near side depends on the interplay between Vi and 
V 2 (higher order terms are negligible). Although the V 2 
parameters are similar, the significantly more negative 
Vi in low- versus high-multiplicity events eliminates the 
near-side V 2 peak in A^. The same applies also to the 
TPC-FTPC correlation comparison between the Au- and 
d-going directions. The V 2 values are rather similar for 
FTPC-d (forward rapidity) and FTPC-Au (backward ra¬ 
pidity) correlations, but the more negative Vi for d-going 
direction eliminates the near-side V 2 peak. If the relevant 
physics in d-|-Au collisions is governed by hydrodynam¬ 
ics, then it may not carry significance whether or not 
there exists a finite near-side long-range correlated yield, 
which would be a simple manifestation of the relative Vi 
and V 2 strengths. 

Our V 2 data are qualitatively consistent with that from 
PHENIX [5D] . While PHENIX focused on the pr depen¬ 
dence, we study the Fourier coefficients as a function of 
A 77 afforded by the large STAR acceptance, as well as the 
event multiplicity. Hydrodynamic effects, if they exist in 
d-|-Au collisions, should naively differ over the measured 
multiplicity range and between Au- and d-going direc¬ 
tions. However, the V 2 parameters are approximately 
constant over multiplicity, and quantitatively similar be¬ 
tween the Au- and d-going directions. On the other 
hand, the correlation comparisons between low- and high- 
activity data reveal different trends for the Au- and d- 
going directions. The high- and low-activity difference in 
the FTPC-Au correlation in Fig.j^b) may resemble ellip¬ 
tic flow, but that in the FTPC-d correlation in Fig. [^c) 
is far from an elliptic flow shape. In combination, these 
data suggest that the finite values of cannot be ex¬ 
clusively explained by hydrodynamic anisotropic flow in 
d-|-Au collisions at RHIC. 

In summary, dihadron angular correlations are re¬ 
ported for d-l-Au collisions at = 200 GeV as a 

function of the event activity from the STAR experi¬ 
ment. The event activity is classified by the measured 
zero-degree neutral energy in ZDC, the charged hadron 
multiplicity in FTPC, both in the Au-going direction, 
or the multiplicity in TPC. In a recent paper we have 
shown that the short-range jet-like correlated yield in¬ 
creases with the event activity [53] ■ In this paper we 
focus on long-range correlations at large lAryj, where jet¬ 
like contributions are minimal on the near side, although 
the away side is still dominated by jet production. Two 
approaches are taken, one to extract the correlated yields 
above a uniform background estimated by the ZYAM 
method, and the other to calculate the Fourier coeffi¬ 
cients, Vn = (cosnAcj)), of the dihadron Acj) correlations. 
The following points are observed: (i) The away-side cor¬ 
related yields are larger in high- than in low-activity col¬ 
lisions in the TPC and FTPC-Au, but lower in FTPC- 
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Uncorrected mid-rapidity dN^^/dri 


FIG. 4: Fourier coefficients (a) Vi, (b) V 2 , and (c) V 3 versus 
the measured mid-rapidity charged particle dNch/dr;. Event 
activity selections by both ZDC-Au and FTPC-Au are shown. 
Trigger particles are from TPC, and associated particles from 
TPC (triangles), FTPC-Au (circles), and FTPC-d (squares), 
respectively. Systematic uncertainties are estimated to be 
10% on Vi and V 2 , and smaller than statistical errors for 
V 3 . Errors shown are the quadratic sum of statistical and 
systematic errors. The dashed curves are to guide the eye. 


d; (ii) Finite near-side correlated yields are observed at 
large Ar/ above the estimated ZYAM background in high- 
activity collisions in both the TPC and FTPC-Au (re¬ 
ferred to as the “ridge”); (iii) The ridge yield appears 
to scale with the away-side correlated yield at the cor¬ 
responding At; < — 1 , which is dominated by the away- 
side jet; (iv) The V 2 coefficient decreases with increasing 
I Ary I, but remains finite at both forward and backward 
rapidities (lAryj « 3) with similar magnitude; (v) The 
Vi coefficient is approximately inversely proportional to 
the event multiplicity, but the V 2 appears to be indepen¬ 
dent of it. While hydrodynamic elliptic flow is not ex¬ 
cluded with a coincidental compensation of jet dilution 
by increasing flow contribution with multiplicity and an 
unexpected equality of elliptic flow between forward and 
backward rapidities, the data suggest that there exists a 
long-range pair-wise correlation in d-l-Au collisions that 
is correlated with dijet production. 


We thank the RHIC Operations Group and RCF at 
BNL, the NERSC Center at LBNL and the Open Sci¬ 
ence Grid consortium for providing resources and sup¬ 
port. This work was supported in part by the Offices 
of NP and HEP within the U.S. DOE Office of Sci¬ 
ence, the U.S. NSF, the Sloan Foundation, the DFG clus¬ 
ter of excellence ‘Origin and Structure of the Universe’ 
of Germany, CNRS/IN2P3, STFC and EPSRC of the 
United Kingdom, FAPESP CNPq of Brazil, Ministry of 
Ed. and Sci. of the Russian Federation, NNSFC, CAS, 
MoST, and MoE of China, GA and MSMT of the Czech 
Republic, FOM and NWO of the Netherlands, DAE, 
DST, and CSIR of India, Polish Ministry of Sci. and 
Higher Ed., Korea Research Foundation, Ministry of Sci., 
Ed. and Sports of the Rep. Of Croatia, Russian Ministry 
of Sci. and Tech, and RosAtom of Russia. 


[1] I. Arsene et al. (BRAHMS Collaboration), Nucl.Phys. 
A757, 1 (2005), nucl-ex/0410020. 

[2] B. Back et al. (PHOBOS Collaboration), Nucl.Phys. 
A757, 28 (2005), nucl-ex/0410022. 

[3] J. Adams et al. (STAR Collaboration), Nucl.Phys. A757, 
102 (2005), nucl-ex/0501009. 

[4] K. Adcox et al. (PHENIX Collaboration), Nucl.Phys. 
A757, 184 (2005), nucl-ex/0410003. 

[5] B. Muller, J. Schukraft, and B. Wyslouch, 

Ann.Rev.Niid.Part.Sci. 62, 361 (2012), 1202.3233. 

[ 6 ] J.-Y. Ollitrault, Phys.Rev. D46, 229 (1992). 

[7] U. Heinz and R. Snellings, Ann.Rev.Nucl.Part.Sci. 63, 
123 (2013), 1301.2826. 

[ 8 ] F. Wang, Prog. Part. Nucl. Phys. 74, 35 (2014), 
1311.4444. 

[9] J. Adams et al. (STAR Collaboration), Phys.Rev.Lett. 
95, 152301 (2005), nucl-ex/0501016. 

[10] B. Abelev et al. (STAR Collaboration), Phys.Rev. C80, 
064912 (2009), 0909.0191. 

[11] B. Alver et al. (PHOBOS Collaboration), Phys.Rev.Lett. 

104, 062301 (2010), 0903.2811. 

[12] B. Abelev et al. (STAR Collaboration), Phys.Rev.Lett. 

105, 022301 (2010), 0912.3977. 

[13] S. Chatrchyan et al. (CMS Collaboration), Phys.Lett. 
B724, 213 (2013), 1305.0609. 

[14] K. Aamodt et al. (ALICE Collaboration), Phys.Rev.Lett. 
107, 032301 (2011), 1105.3865. 

[15] B. Alver and G. Roland, Phys.Rev. C81, 054905 (2010), 
erratum-ibid. C82, 039903 (2010), 1003.0194. 

[16] L. Adamczyk et al. (STAR Collaboration), Phys. Rev. 
C 88 , 014904 (2013), 1301.2187. 

[17] V. Khachatryan et al. (CMS Collaboration), JHEP 09, 
091 (2010), 1009.4122. 

[18] S. Chatrchyan et al. (CMS Collaboration), Phys.Lett. 
B718, 795 (2013), 1210.5482. 

[19] B. Abelev et al. (ALICE Collaboration), Phys.Lett. 
B719, 29 (2013), 1212.2001. 

[20] G. Aad et al. (ATLAS Collaboration), Phys.Rev.Lett. 
110, 182302 (2013), 1212.5198. 

[21] B. B. Abelev et al. (ALICE Collaboration), Phys.Lett. 
B726, 164 (2013), 1307.3237. 










[22] P. Bozek, Eur.Phys.J. C71, 1530 (2011), 1010.0405. 

[23] P. Bozek and W. Broniowski, Phys.Lett. B718, 1557 
(2013), 1211.0845. 

[24] A. Dumitm et ak, Phys.Lett. B697, 21 (2011), 

1009.5295. 

[25] K. Dusling and R. Venugopalan, Phys.Rev. D87, 054014 
(2013), 1211.3701. 

[26] K. Dusling and R. Venugopalan, Phys.Rev. D87, 094034 
(2013), 1302.7018. 

[27] D. Molnar, F. Wang, and C. H. Greene (2014), 1404.4119. 

[28] A. Adare et al. (PHENIX Collaboration), Phys.Rev.Lett. 
Ill, 212301 (2013), 1303.1794. 

[29] L. Adamczyk et al. (STAR Collaboration) (2014), ac¬ 
cepted by Phys.Lett.B, 1412.8437. 

[30] A. Adare et al. (PHENIX Collaboration) (2014), 
1404.7461. 

[31] J. Adams et al. (STAR Collaboration), Phys.Rev.Lett. 
91, 072304 (2003), nucl-ex/0306024. 

[32] B. Abelev et al. (STAR Collaboration), Phys.Rev. C79, 


034909 (2009), 0808.2041. 

[33] K. Ackermann et al. (STAR Collaboration), 
Nucl.Instrum.Meth. A499, 624 (2003). 

[34] C. Adler et ah, Nucl.lustrum.Meth. A499, 433 (2003). 

[35] M. Anderson et ah, Nucl.Instrum.Meth. A499, 659 
(2003), nucl-ex/0301015. 

[36] K. Ackermann et ah, Nucl.Instrum.Meth. A499, 713 
(2003), nucl-ex/0211014. 

[37] N. Ajitanand et ah, Phys.Rev. C72, 011902 (2005), nucl- 
ex/0501025. 

[38] M. Gyulassy and X.-N. Wang, Comput.Phys.Commun. 
83, 307 (1994), nucl-th/9502021. 

[39] L. Xu, C.-H. Chen, and F. Wang, Phys.Rev. C88, 064907 
(2013), 1304.8120. 

[40] G. Aad et al. (ATLAS Collaboration), Phys.Lett. B725, 
60 (2013), 1303.2084. 

[41] B. B. Abelev et al. (ALICE Collaboration), Phys.Rev. 
C90, 054901 (2014), 1406.2474. 



